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THE  THEORY  OF  ROTATIONALLY 
SYMMETRIC  PLASTIC  SHELLS 


By  Philip  G.  Hodge,  Jr. 


ABSTRACT 

The  defining  equations  for  a  rigid/perfectly-plastic  shell 
are  derived  from  basic  principles.  On  the  basis  of  a  single  geo¬ 
metric  assumption  for  the  velocity  field,  generalized  strain  rates 
and  stresses  are  defined  and  equilibrium  relations  deduced.  Shell 
yield  conditions  and  the  flow  law  are  discussed  in  general  terms  and 
then  specifically  for  piecewise  linear  yield  conditions. 

Preceding  the  general  shell  problem*  the  theory  of  beams 
under  bending  and  axial  forces  is  discussed  to  give  a  general  in¬ 
sight  into  plastic  structural  behavior.  The  paper  closes  with  an 
application  to  cylindrical  shells  and  a  discussion  of  areas  for  future 
development. 


1.  INTRODUCTION 


For  historical  reasons ,  the  structural  theory  of  shells  has,  until 
recently,  been  exclusively  an  elastic  theory.  Thus,  in  the  usual  develop¬ 
ment,  the  assumptions  of  a  linear  theory  of  elasticity  and  of  shell  theory 
have  been  introduced  at  the  outset,  and  equations  have  been  derived  which 
apparently  depend  upon  both  of  these  sets  of  assumptions.  As  a  result, 
it  is  desirable  to  begin  an  inelastic  theory  of  shells  by  a  consideration  of 
basic  principles.  In  the  interests  of  possible  other  inelastic  theories,  the 
basic  kinematic  and  static  equations  of  shell  theory  have  been  derived  in 
Sec.  3  by  a  method  which  is  independent  of  any  particular  constitutive 
equations  of  the  material.  The  final  results  will,  of  course,  come  as 
no  surprise  to  those  familar  with  classical  elastic  shell  theory,  but  it  is 
not  obvious  that  the  full  extent  of  their  generality  has  always  been  appreci¬ 
ated  previously. 

Although  aiming  for  maximum  generality  of  material  assumptions, 
it  has  seemed  desirable  to  limit  the  treatment  here  to  a  rotationally  symmetri 
shell  under  rotationally  symmetric  loading  with  no  torque.  This  limitation 
considerably  simplifies  the  mathematical  presentation  of  the  fundamental 
concepts  presented,  although  the  concepts  themselves  are  evidently  appli¬ 
cable  to  the  more  general  case  of  arbitrary  shells  under  arbitrary  loadings. 

Further  justification  lies  in  the  fact  that  very  few  problems  have  been  solved 

$ 

to  date  which  do  not  fall  into  the  rotationally  symmetric  catagory.  For 

*  The  only  two  exceptions  known  to  the  author  are  a  paper  by  Fialkow  [  1] 
in  which  bounds  are  found  on  a  cylindrical  roof,  and  one  by  Hodge  and 
Panarelli  [  2]  concerned  with  a  cylindrical  shell  under  pressure,  end  load, 
and  torque.  (Numbers  in  brackets  refer  to  the  list  of  references  collected  at 
the  end  of  the  paper). 
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brevity  of  notation,  we  shall  use  the  unmodified  term  "shell"  to  refer  to 
the  fully  rotationally  symmetric  catagory  defined  above* 

By  way  of  introduction  to  the  shell  problem,  we  present  in 
Sec.  2  a  development  of  the  general  and  plastic  theory  of  beams  under 
combined  bending  and  axial  force  which  was  first  formulated  by  Onat 
and  Prager  [3,  4]  .  In  this  much  simpler  case,  we  can  show  in  some 
detail  how  a  structural  theory  is  derived  from  a  three-dimensional  theory. 
A  single  kinematical  assumption  is  made  concerning  the  deformation  or 
velocity  field.  Using  this  assumption  and  the  principle  of  virtual  work, 
one  is  led  naturally  to  the  definition  of  generalized  strain  rates  and 
generalized  stresses.  As  first  pointed  out  by  Pf  agfef  f  5}  ,  thethecrrernB 
of  limit  analysis  apply  immediately  to  such  variables;  Hodge  [  6,  7]  has 
shown  that  the  general  theory  of  elasticity  and  various  models  of  plasti¬ 
city  can  all  be  conveniently  developed  in  terms  of  generalized  variables. 

When  the  constitutive  equations  of  a  rigid/perfectly-plastic 
material  are  inserted  into  the  beam  problem,  the  concept  and  principal 
properties  of  the  yield  curve  and  plastic  flow  law  are  all  derived  in  a 
simple  manner.  That  these  properties  are  all  valid  for  any  number'  of 
generalized  strain-rate  and  stress  variables  is  certainly  plausible.  Al¬ 
though  they  can  all  be  proved  based  upon  Drucker's  "stability  postulate" 
[8,  9]  ,  it  has  seemed  more  in  keeping  with  the  spirit  of  the  present  paper 
to  assume  them  and  develop  from  them  the  specific  plastic  constitutive 
equations  for  the  shell  problem,  as  is  done  in  Sec.  4. 

The  derivation  of  yield  conditions  for  plastic  shells  has 

$ 

been  the  subject  of  numerous  papers.  Although,  the  resulting  yield 
conditions  differ  greatly  in  mathematical  complexity,  the  theorems  of 


4  An  account  of  most  of  these,  together  with  original  source  references, 
can  be  found  in  [  10]  . 
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limit  analysis  can  be  invoked  to  bound  their  differences  in  predictions* 
Since  very  little  direct  experimental  evidence  is  available  upon  which 
to  base  a  reasonable  choice*  it  has  seemed  advisable  to  concentrate 
attention  upon  the  mathematically  simple  class  of  piecewise -linear 
yield  conditions  in  Sec.  5.  Not  only  does  this  greatly  simplify  the 
solution  of  problems*  but*  as  is  shown*  the  results  can  be  used  to 
bound  the  yield-point  loads  of  any  conceivable  material  yield  condition. 

The  primary  purpose  of  the  present  paper  is  to  survey  the 
theory  of  plastic  shells*  rather  than  to  give  a  catalogue  of  problem 
solutions.  However*  in  the  interest  of  providing  a  typical  application 
which  is  mathematically  simple*  we  have  included  an  example  concerned 
with  a  cylindrical  shell  in  Sec.  6. 

The  treatment  in  Sections  4  through  6  was  intentionally  res¬ 
tricted  to  an  idealized  rigid/perfectly-plastic  material.  The  final  section 
of  the  paper  discusses  the  reasons  for  thiB  limitation  and  indicates  quali¬ 
tatively  some  possible  generalizations. 


2.  BEAM  THEORY 


As  a  prototype  of  shell  theory,  we  shall  first  examine  the 
much  simpler  case  of  an  initially  straight  beam  under  the  influence  of 
bending  moment  M  and  direct  stress  N.  We  shall  make  the  initial 
assumption  that  plane  normal  sections  remain  inextensible,  plane,  and 
normal  to  the  centroidal  aris.  At  any  cross  section  x,  then  the  velo¬ 
city  field  is 

u  =  V  -  z  w'  u  =  0  u  =  W  (2.1) 

x  y  z 

wheje  V  and  W  are  functions  of  x« 


The  strain-rate  tensor  field  associated  with  (2.  1)  has  the 
axial  strain  as  its  only  non-vanishing  component: 

€  =  V  -  zw"  (2.2) 

x 


so  that  the  total  internal  rate  of  work  reduces  to 

wint  -  lv’u  Uj  dv  - 1 1  ** (v'  - 1  w"»  dA  d* 


(2.3) 


Since  V  and  W  depend  only  on  x,  we  may  rewrite  (2.  3)  in  the  form 


W.  .  =  /  <NV'  -  MW ")  dx 
int  0 


(2.4) 


where 


N  =  /  <r  dA 
A  X 


M  = /  a  zd A 
A  x 


(2.5) 


Since  N  and  M  are  the  stress  resultants  for  direct  stress 
and  bending  moments,  respectively,  they  are  an  obvious  choice  for 
generalized  stresses  for  the  beam  problem.  It  then  follows  from  Prager 
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[5]  that  the  generalized  strain  rate  should  be 

■  ii 

e  =  V  K  =-W 

whence  we  can  rewrite  (2.  2)  in  the  form 

f.  =  e  +  z  K 
x 


(2.6) 


(2.7) 


An  alternative  form  of  (2.  4)  is  obtained  by  integrating  it 
twice  by  parts: 

1  i  n  ■  i  l-* 

(2.8) 


W.  =  -  /  (n'  V  +  m"w)  dx  +  [NV  +  MW  -  MW] 


int 


The  external  rate  of  work  may  be  written 

L  ,  L 

W  =  /  PWdx  +  [5lV  +  SW  -  MW  ] 

©Xt  A  A 


(2.9) 


where  P  is  the  applied  normal  load  and  N,  S,  and  M  are,  respectively, 
the  axial  force,  shear  force,  and  bending  moments  applied  at  the  beam 
ends. 


The  principle  of  virtual  work  rate  states  that  the  internal  and 
external  work  rates  must  be  equal  for  all  sufficiently  continuous  velocity 
fields  V,  W.  Obviously  this  condition  requires 

N  =0  m"  +  P  =  0  0<  x<L 


(2.  10) 

N  =  N  M=S  M  =  M  x  =  0,  x  =  L 


The  development  thus  far  is  independent  of  any  material  property. 
Since  only  the  axial  force  <r  enters  in  the  definition  of  the  generalized 
stresses,  we  shall  assume  that  it  is  the  only  stress  to  influence  the  mater¬ 
ial  behavior.  In  particular,  for  a  rigid/ perfectly-plastic  material,  the 
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the  yield  behavior  is  fully  characterized  by  the  yield  stress  o^. 

If  =  <Tq  at  an  element,  the  element  may  have  any  positive  axial 
strain  rate  e^;  if  <7^  =  -  it  may  have  any  negative  axial  strain  rate 
rate;  if  *°o<<Tx  <O0  t*le  e^ement  must  remain  rigid;  stress  states 
lxr0  are  not  tolerated.  In  other  words  the  strain  rate  can  be 
non-zero  only  if  <7  =  Oq  for  positive  and  <7  =  -  for  negative  €*. 

Since  (2.  7)  describes  a  strain  rate  which  is  linear  in  z,  there 
will  be  a  particular  value  z  *  £  H  at  which  =  0.  For  all  z  on  one  side 

£  H,  £  will  be  positive  and  hence  <7  =  trn,  whereas  on  the  other 

XU 

side«  **  w*W  be  ne£atiye  an(J  h$»ce  <^=  -  fy.  ABBUmift&  fox  dcfinite  - 
ness  that  K  is  positive,  we  see  that  the  strain-rate  and  stress 
distributions  must  have  one  of  the  forms  shown  in  Fig.  1.  Therefore, 
it  follows  from  (2.5),  (2.  6),  and  Fig.  1  that 

M  =  0 

£<-  1  (2.  11a) 

N  a  <7  [  -  /  B(z)  dz  +  /  B(z)  dz] 

-H  £  H 

£  H  H 

M  =<7.  [  -/  z  B(z)dz  +  /  z  B  (z)  dz] 

-H  £  H 

e/( -  £H)  =  K>0  for  -  ls£s=l  (2.11b) 

H 

N  =  -  Oq  /  B  (z)  dz  =  -  l^j  M  =  0 
-H 

-e/HSKaO  for  Is  £  (2.11c) 


H 

N  =  (7^  /  B(z)  dz  =  I^J 
-H 

e/H>K>0  for 

£H  H  . 


where  B(z)  is  the  width  of  the  section 
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Equation  (2.11b)  defines  a  curve  in  the  N,  M  plane.  If  we  consider 
the  case  K<  0,  •  1  s  J  <  1,  we  obtain  the  image  in  the  origin  of  (2.11b). 
Evidently  the  sum  of  these  two  curves  is  closed.  This  resulting  closed 
curve  is  called  the  yield  curve.  If  the  stress  point  N .  M  is  on  the  yield 
curvef  plastic  flow  can  take  place;  if  the  stress  point  is  inside  the  curve, 
the  section  is  rigid;  stress  points  outside  of  the  yield  curve  are  not  tolorated. 
Figure  2  shows  the  resulting  yield  curve  for  a  rectangular  section. 

It  follows  from  (2.11b)  that 

dM/dN  =  £  H  (2.12) 

whence  the  yield  curve  is  evidently  convex.  Further*  since 

K/e  =  -  \U  H  (2.  13) 

the  "strain-rate  vector"  with  components  (e.  K)  is  normal  to  the  yield 
curve  at  the  stress  point  i  . 

Equations  (2.11a)  give  the  single  stress  point  (I^»  0)  but  permit 
a  variety  of  strain-rate  vectors.  If  i  tends  to  -1  with  positive  K,  then 
it  follows  from  (2.  13)  that 

lim  (e,  K)  =  X.  (1,  1/H)  (2.  14a) 

i  1  1 

where  X^  is  an  arbitrary  positive  scalar.  Similarly,  if  K  is  negative. 

lim  (e,  K)  *  X,.  (1,  -  1/H)  (2.14b) 

t  +  - 1  C 

where  X^  a 0.  From  the  inequality  (2.  11a)  on  e/H  and  K.  together  with  the 
corresponding  inequality 


e/H  a  -  K  a  0 
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for  negative  K,  we  see  that  at  £  =  -  1,  e  and  K  are  restricted  only 
by 

e/H  >  IK  |s  0  (2.15) 

Finally,  it  is  obvious  that  the  sum  of  the  two  vectors  (2.14)  will  satisfy 
(2.15)  for  any  non-negative  choice  of  Xj  and  X^.  Typical  strain-rate 
vectors  are  shown  in  Fig.  2  for  a  regular  point  B  governed  by  (2.11b) 
and  for  a  singular  point  A  governed  by  (2.  11a). 

For  the  beam  problem  considered  here,  we  have  started  with  a 
simple  geometrical  assumption  and  with  the  tensile  behavior  of  the  material. 
We  have  then  developed  the  equilibrium  conditions,  the  convexity  of  the 
yield  curve,  the  normality  of  the  strain-rate  vector  of  a  regular  point, 
and  the  behavior  of  the  strain-rate  vector  at  a  singular  point.  It  can  be 
shown  that  these  concepts  are  all  carried  over  to  the  more  general  problem 
of  rotationally  symmetric  plastic  shells.  We  shall  use  these  facts  in  the 
next  two  sections. 


3.  BASIC  EQUATIONS  OF  SHELL  THEORY 


The  theory  of  elastic  shells  has  been  thoroughly  studied,  whereas 
investigations  of  inelastic  shells  are  quite  recent.  As  a  result,  it  is 
not  always  clear  which  aspects  of  elastic  shell  theory  are  dependent  upon 
the  elastic  constitutive  equations  and  which  are  equally  valid  for  other 
materials.  Therefore,  we  shall  begin  this  discussion  by  deriving  the 
fundamental  rotationally  symmetric  shell  equations  in  a  manner  which  makes 
no  reference  to  any  constitutive  equations.  The  results  will  thus  be  appli¬ 
cable  to  shells  of  any  material  and  hence,  in  particular,  to  a  rigid/perfectly- 
plastic  material. 

A  single  kinematic  assumption  is  made  concerning  the  type  of  de¬ 
formation  which  the  shell  undergoes,  namely: 

Initially  straight  normals  to 

the  undeformed  middle  surface 

remain  straight,  inextensible,  (3.  1) 

and  normal  to  the  deformed 

middle  surface. 

The  degree  of  validity  of  this  assumption  as  applied  to  a  real  physical  struc¬ 
ture  determines  the  validity  of  calling  the.  structure  a  shell  in  the  sense  used 
in  this  paper.  The  validity  of  (3.  1)  cannot  be  checked  internally  on  the  basis 
of  shell  theory  but  must  be  done  on  the  basis  of  complete  three -dimensioned 
solutions,  experimental  evidence,  and/or  sound  engineering  judgement.  It  is 
certainly  conceivable  that  its  validity  should  depend  not  only  on  the  geometrical 
parameters  such  as  thickness  and  radii  of  curvature,  but  on  the  type  of  mater¬ 
ial  being  considered,  i.  e. ,  elastic,  plastic,  etc. 
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A  detailed  discussion  of  the  situations  under  which  (3.  1)  is  a 
reasonable  approximation  to  reality  would  be  beyond  the  scope  of  the 
present  paper,  even  if  it  could  be  done  in  any  generality.  Therefore,  in 
the  remainder  of  the  present  development  we  shall  assume  (3.  1)  and  in¬ 
vestigate  its  consequences. 

The  method  of  attack  is  similar  to  that  used  in  the  previous  section. 
We  shall  first  find  the  expressions  for  linear  strain-rate  components 
for  a  general  deformation  of  a  shell.  Next  we  particularize  these  strain- 
rate  velocity  relations  in  accord  with  assumption  (3.  1).  The  internal 
work  rate  associated  with  an  arbitrary  displacement  field  subject  to  (3.1) 
is  then  computed  and  shown  to  naturally  suggest  the  generalized  strain- 
rate  variables  appropriate  to  the  rotationally  symmetric  shell  problem. 
Finally,  the  principle  of  virtual  work  is  invoked  to  define  the  generalized 
stress  variables  and  to  obtain  the  equilibrium  equations  which  they  must 
satisfy. 

We  consider  an  arbitrary  longitudinal  plane  of  the  shell  and  let  O 
be  an  arbitrary  point  on  the  center  line.  Fig.  3.  The  coordinate  <l>  is 
defined  as  the  angle  between  the  axis  and  a  normal  to  the  center  line  through 
O;  the  coordinate  4  is  measured  inward  along  the  normal  from  the  center 
line.  Let  points  P  and  R  have  the  same  coordinate  as  O  with  4 
coordinates  4  and  4  +  A4  ,  respectively;  let  Q  have  coordinates 
<p  +  A  <p  and  4  ,  all  as  shown  in  Fig.  3. 

It  is  now  convenient  to  introduce  a  set  of  Cartesian  coordinates 
r,  z  in  the  plane  of  Fig.  3.  The  following  vectors  may  be  readily  identified: 
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PQ  *  (Rj  -  {  )A  <t>  (cos  <t>,  sin  <p) 
PR  aA|  (-sin  <p,  cos  0) 


(3.2) 


Consider  next  a  small  vector  velocity  field  .  If  and 
are  the  components  of  ^u,  in  the  £  and  <t>  directions,  then  maybe  re¬ 
ferred  to  r  and  z  components  by 

ju  =  un  (-sin  <p  ,  cob  <t>  )  +  u^  (cos  <t>,  sin  <t>)  (3.3) 

i  i  i 

Let  P|  Q,  and  R  denote  the  respective  positions  of  P,  Q,  and  R  after 
the  displacement.  Then 


P  Q  =  PQ  +  QQ  -  PP  =  PQ  +  uJQ)  -  u  (P) 
hence,  in  view  of  Eqs.  (3.  2)  and  (3.  3) 

P  Q  =  (Rj  -  £  )  A  <p[(l  +  ^ . n  )  (cos  <f>,  sin  <p)  +  (u^  +  uj  )  (-sintf>,  cos  <p)] 


(3.4) 


PR=A{[(I+u)  (-sin  <p,  cos  <p)  +  u.  (cos  <f>,  sin  <t>  )] 
n  v* 


Here  we  have  regarded  A  <f>  and  A  4  as  infinitesimals  and  have  denoted 
differentiation  with  respect  to  <p  and  |  by  a  prime  and  a  dot,  respectively. 

If  l^u/Rjl,  Iju/Rjl  ,  and  IJi  I  are  all  small  compared  to  unity  we  can 
easily  deduce  from  (3.4)  and  (3.  2)  that 


i  i 


IPQI  =  (R1  -  £  )*<t>(l+ 


PR  \  =  A  £  (1  +  u,  ) 

n 


(3.5) 


f  t 


P  Q  PR®  A<#>A£  (Rj  -4  )  (h0+ 


PQ  I  =  (Rj  -  4  )A 


I  PR  I  =  A-( 


PQ  •  PR  =  0 
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Extensional  strain  rate  is  defined  as  the  rate  of  change  in  length 
per  unit  length  of  a  line  element.  Therefore  the  extensional  strain  rate 
of  an  element  in  the  direction  of  p£  is 


P  R 


PR  I 


n 


and,  in  similar  fashion 


'■<f> 


I  PR  I 


PQ1  I  -  I  PQ I 


=  u 


"POT 


n 


=  j£_  Un 


(3.  6a) 


VT 


(3.  6b) 


The  extensional  strain  rare  in  th  e  cireamrer  eiiliai  air  ectlfih  is  rare  FaTe  aT 
increase  in  the  circumference  of  a  circle  through  P  divided  by  the  original 
circumference.  Since  the  circumference  is  proportional  to  the  radius,  we 
may  equally  well  use  the  r  components  of  the  vectors  CP  and  CP  to 
define 


€*  = 


(CP)r  -  <CP)r 


<C?) 


=  U«PCOB0  ~Un9in^(3.6c) 
(Rg  -  {  )  sin  0 


The  shear  strain  rate  is  defined  as  the  rate  of  change  in  an 
initially  right  angle.  In  view  of  the  restriction  to  rotational  symmetry 

=r*i  ■»  <3-M> 


whereas  the  remaining  shear  strain  rate  is 


ii  ii 

„  _  P  Q  *  P  R 

n<P  I  PQ1!  I  p'r  I 


=  u 


UA  ,  u  ' 

I  *  ±  n 


(3.6e) 
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Equations  (3.  6)  give  the  strain  rate  components  at  any  point  in 
the  shell  for  an  arbitrary  displacement  field.  We  now  introduce  the 
assumption  (3.  1).  If  Vn  and  represent  the  velocities  of  a  point  O 
on  the  center  line,  then  it  follows  from  Fig.  4  that  at  P 

”„-v„  \  <3-7' 
Substituting  (3.  7)  into  (3.  6)  we  find  that  for  the  particular  velocity 
field  being  considered, 

€n  a  7n*  *  =0 


(3.  8) 


The  internal  rate  of  work  associated  with  a  rotationally  symmetric 

shell  may  be  written 

_  ot  J 

wi».  ■  a  § 

0  -H  0 

where  the  shell  extends  from  <p  =  0  to  <f>  =<*  and  is  of  thickness  2H;  here 


c2r 

J  wint  ^R1  "  ^  sin  <pd9d$d<p  (3.9) 


wint  =  ffn  €n  +  <r0  €0  +  ^  *  <f> +Tn<l>yn<t>+ T 9  ^9  *  +TJnT«n 

We  substitute  the  particular  strain  rate  field  (3. 8)  into  (3.  9)  arid  note  that 
V  V  Rj.and  are  independent  of  £  and  9  to  write  the  resulting  ex¬ 
pression  in  the  form 
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[' 5"  -H  i  4 


«2 


*  f<2 


sin  <£d  <f> 


(3.  10) 


The  factors  in  the  first  bracket  of  each  term  of  {3.  10)  are  chosen  as 
our  generalized  strain  rates.  They  may  be  given  simple  physical  interpreta¬ 
tions  as  extension  rates  and  approximate  curvature  rates  of  the  middle  surface. 
Thus 


•*  *  <v*  «*  ♦  -  V  ih.  v“(V  '  V»,/RI 


(3.11a) 

As  first  pointed  out  by  Prager  [5]  >  generalized  stresses  and  strain 
rates  must  be  chosen  so  that  the  internal  work  rate  is  proportional  to  their 
scalar  product.  Therefore*  we  define  the  second  bracketed  factor  of  each 
term  as  a  generalized  stress: 
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(3.  lib) 


where  we  have  assumed  for  simplicity  that  the  stress  distribution  also 
is  rotationally  symmetric. 

We  observe  that  (3.  11)  are  precisely  the  total  force  per  unit  length 
and  total  moment  per  unit  length  about  the  center  line  which  are  trans- 
rhltfed  across  an  element  ofihe  shell.  Therefore,  the  definitions  T3.  ITT 
are  in  agreement  with  the  usual  definitions  of  stress  resultants. 


External  work  may  be  done  on  the  shell  by  surface  loads  applied 
at  {  *  +  H,  by  edge  loads  applied  at  <t>  =  P  and  <l>  =«t ,  or  by  body 
forces.  We  denote  by  Tn  and  the  force  per  unit  area  components  trans¬ 
mitted  across  a  surface  4  *  const,  from  greater  to  lesser  values  of  4  » 
by  Sn  and  similar  forces  across  a  surface  <l>  =  const.,  and  by  Fn  and 
the  body  force  per  unit  volume  components.  Then  the  external  rate  at 
which  work  is  done  during  an  arbitrary  symmetric  deformation  is 


2ir  H 


Wext  I  M  <Sn"n+S*V 

0  -H 


* 

P 


(3.  12) 


+  t  J  <Tnun  +  T*V  (Rl“fi)  (R2 -«)  sin«d<#»]H 

P 


-H 


t  -H  (Fnun  +  )  (Rj  -  ( )  (Rj  -{, )  .in  *d{  d* 


mssmmmm 


afai  ii  m&m  &bi>  a.  ■«;  a&a  asfet  >. .  &fc^ 
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We  substitute  the  particular  velocity  field  {3.  7)  into  (3.  12)  and  re¬ 


arrange  the  terms  to  obtain 

>  .rf,  /f  s  h—i  d()  +  H 

f-H  "  *2  '  «.h  ' 


v.  +  v  •  /.F 

VYi“ 


0  R- 


jj  *—]- 


+{F-^ 


-  «  *2  > 

Vdv 


+  V 
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(fvV 


-  H 

H 

4  7 

H 

-J 

F0 

-H 

-H 

H  H 

V* 

7  +)  ‘ 

«Z  . 

L  -H 

Ri-I  R2  -4 


d0 


VV/r  Ri_l  v«7  +  C  aF  ri-«  «»-« \) 

— —  ( L€T^ -^r — *rJH  i  *  *  “iq — sin^d 


(3.  13) 

Just  as  the  factors  in  (3.  10)  were  recognized  as  resultant  internal  forces 
and  moments,  so  the  factors  in  parentheses  in  (3.  13)  may  be  recognized 
as  resultant  applied  forces  and  couples.  Thus  we  are  led  to  define  bound¬ 
ary  forces 


=  !Sn 


•H  *2-« 
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■d« 


R, 
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-H 
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■I 

-H 


*2  “  * 

S0  -V"d* 
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=  |  d| 


(3.  14) 


0  "  J  s"0  ~RT 

-H  £ 


at  the  edges  0  =  (3  and  0  =o<  of  the  shell,  and  distributed  loads  and 
couple  by 
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(3.  14b) 


for  0  <  <t><*  .  The  definitions  (3.  14)  are  consistent  with  the  usual 
physically  based  definitions  and  include  the  possibility  of  a  couple 
being  distributed  over  the  surface. 


We  next  set  the  internal  rate  of  work  (3.  10)  equal  to  the  external 
rate  of  work  (3.  13),  using  the  definitions  (3.  11)  and  (3.  14).  Integrat¬ 
ing  the  term  involving  by  parts  we  can  write  the  resulting  equation 
in  the  form 

zr  <we«  -  wmt' '  J  {V-R1  N»  “•♦  +  Rl,*z  V*”** 

0 


+  Vn  (Rj  Ng  sin  «#>  +  Rg  sin  ^  +  Rj  Rg  Pn  sin  <p ) 


V±  +  V_ 


-  p  Rj  Mq  cob  4>  +  (Rg  sin  <t>)  +  Rj  Rg  sin  <*>]jd  <f> 

+  t<v»s  +  v*V  ~ *  ft,---  <*V*V}  R2*ia*|' 


(3.  15) 


18 


The  group  of  terms  in  the  bracket  in  (3.  15)  suggest  the  definition  of  a 
quantity  S  by 

Rj  R2  S  sin#  s^^M^sin#)  -  Rj  cos  <f>  +  sin#  (3.  16) 

whence  a  further  integration  by  parts  of  (3.  15)  leads  to 

j*  ^[<*2  N#  Bin<t>)'  ‘  R1  N0  sin#  +  Rj  P^sin #] 

P 

+  Vn  [(R^  S  sin#)'  +  Rj  sin#  +  Rg  sin#+.Rj  R2  Pn  sin#)}  d# 

v+v.  *  <3‘17> 
*  [  (Vn  (S  -  5)  +  V^Oylty  -  It,  sin#]^ 

If  Eq.  (3.17)  is  to  be  valid  for  all  sufficiently  continuous  velocity  functions 
and  V^»  then  the  two  factors  in  square  brackets  on  the  left-hand  side 
must  vanish  indentically  in  #  and  the  three  parenthentical  factors  on  the 
right-hand  side  must  vanish  at  the  shell  edges.  Thus 


(1*2  N(f>  Bin  ♦)'  -  Rj  cos  #  -  Rg  S  sin  #  +  Rj  R2  P^  sin  #  =  0 

(R2  S  sin  #)'  +  Rj  Ng  sin  #  +  R 2  sin  #  +  Rj  P^  sin  #  »  0 
for  all  #(  and  at  #  =ei  and  #  =  p 


(3.18) 


8  =  5  N#  =  N# 


(3.  19) 


Equations  (3.  16)  and  (3.  18)  are  the  familiar  equations  of  equilibrium 
for  a  rotationally  symmetric  shell  (in  the  usual  case  =  0),  and  (3.19) 
express  the  fact  that  shear  force,  normal  direct  stress,  and  normal 
bending  moment  must  each  be  continuous  at  the  boundary. 
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In  the  following,  it  will  prove  convenient  to  deal  exclusively 

with  dimensionless  quantities.  To  this  end  we  denote  by  Nq  the 
maximum  direct  stress  which  the  shell  can  withstand  in  uniaxial  tension 
and  by  Mq  the  maximum  uniaxial  bending  moment  and  define 

n*  s  N*  '  No  s  N9 '  *0 

™OaMetM0  (3.20a) 

Ke  -oyvt  K9 

Further  we  let  A  represent  a  typical  dimension  of  the  shell  and  define 


r  =  R/A 

r,  •  Rj/A 

r2  =  Pg/A 

v  =  V  /A 
n  n 

>"VA 

h  =Mo/a^ 

Pn  =  PnA/N0 

P,  *  P<,  A'^ 

s  =S/1^ 

'V  =  wta./2'V2 

%X.“W.«/2'^A2 

In  terms  of  these  dimensionless  quantities,  the  generalized  strain 
rates  and  stresses  (3.  11)  become 
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and  the  equilibrium  equations  (3*  16)  and  (3.  18)  (with  *  0)  may  be 
written 

"  rln0  cos  *  ‘  r8  +  rt\P<f>  e  0 
(rs)'  +  r^n^  sin  <t>  +  rn^  +  rriPn  =  ®  (3.22) 

h  [  (rir^)1  -  r^nig  cos  <t>]  -  rr^e  =  0 


where 


r  =  r^  sin  0 


(3.23) 


is  the  dimensionless  distance  from  the  axis. 


4.  CONSTITUTIVE  EQUATIONS  OF 
PLASTIC  SHELL  THEORY 


The  basic  unknowns  of  shell  theory  are  the  four  generalized 
stresses  m^,  the  shear  stress  s,  and  the  velocity 

components  and  vn>  Since  Eqs.  (3.  16)  and  (3.  18)  provide  only 
three  equations  for  the  seven  unknowns,  it  is  necessary  to  provide 
four  more  equations  in  these  same  unknowns  which  characterize 
the  particular  shell  material. 

For  a  linear  elastic  material  these  constitutive  equations  are 
trivally  derived.  Hooke's  law  is  substituted  into  (3.  lib)  to  give  general¬ 
ized  stresses  in  terms  of  material  strains,  and  these  latter  are  expressed 
in  terms  of  displacements  by  expressions  analogous  to  (3.8).  It  is  time¬ 
saving  to  express  the  results  in  terms  of  generalized  strains  by  utiliz¬ 
ing  the  displacement  form  of  (3.  21a). 

For  the  plastic  material,  it  is  necessary  to  first  express  the 
yield  condition  in  termB  of  the  generalized  stresses.  For  any  assumed 
form  of  the  material  yield  condition,  the  shell  yield  condition  can  be 
found,  at  least  in  theory.  Assuming  that  H/Rj  and  H/R2  are  negligible 
compared  to  unity,  Onat  and  Prager  [  1 1]  have  found  the  yield  condition 
for  a  uniform  shell  whose  material  satisfies  the  Tresca  yield  condition. 
Similar  treatments  of  the  Mises  yield  condition  applied  to  shell  problems 
have  been  carried  out  by  Hodge  and  Panarelli  [  12,  2]  .  However,  the 
results  are  so  complex  as  to  have  been  relatively  little  used  in  the  solution 
of  problems.  Among  the  few  complete  solutions  (as  opposed  to  bounds 
found  by  limit  analysis)  are  those  by  Hopkins  and  Wang  [  13]  for  a 
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flat  plate,  Hodge  and  Sawczuk  [  14]  for  a  cylindrical  shell,  and  Onat 
and  Lance  [  15,  16]  for  a  nearly  flat  cone. 

In  view  of  the  difficulty  of  obtaining  complete  solutions  ac¬ 
cording  to  the  Tresca  or  Mises  uniform  shell  yield  conditions, 
considerable  interest  has  been  expressed  in  the  use  of  approximate 
yield  conditions.  An  approximate  yield  condition  can  be  viewed  in 
either  of  two  lights.  On  the  one  hand,  by  means  of  the  Bounding 
Surface  Lemma  [  10]  of  limit  analysis,  it  may  be  used  to  provide 
upper  and  lower  bounds  on  the  yield-point  load  according  to  some 
mom  exactyicld  condition.  r,An».rai  factor* 

found  by  Hodge  and  Sankaranarayanan  [  17,  18]  .  Alternatively,  one  may 

visualize  an  ideal  shell  made  of  a  modified  material  such  that  the  approxi 
mate  yield  condition  for  the  real  shell  becomes  the  exact  yield  condition 
for  the  ideal  shell.  If  the  material  difference  between  the  real  and  ideal 
shell  are  not  too  great,  then  the  behavior  of  the  ideal  shell  should  pro¬ 
vide  valuable  information  concerning  the  behavior  of  the  real  shell. 

Further  justification  for  the  use  of  approximate  yield  conditions 
for  shells  is  derived  from  the  fact  that  either  the  Tresca  or  Mises 
material  yield  condition  already  represents  an  approximation  to  reality, 
and  hence  it  is  misleading  to  speak  of  a  shell  yield  condition  based  on 
either  of  them  as  "e-xact". 

Once  the  yield  condition  has  been  decided  upon,  it  provides  one 
of  the  necessary  four  constitutive  equations.  The  remaining  equations 
are  obtained  from  the  flow  law.  The  plastic  potential  flow  law  was 
first  proposed  by  Mises  [  19]  and  later  shown  to  be  a  consequence 
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of  Drucker's  postulates  for  a  stable  material  [8,  9]  .  In  terms  of 
the  generalized  stresses  and  strains  defined  by  (3.21),  we  may 
represent  the  yield  condition  as  a  surface  in  a  four -dimensional 
generalized  stress  space  in  the  form 

f  (nfl  »  m0  ,  =  1  (4.  1) 

The  flow  law  then  states  that  the  strain  rate-vector 

J3  =  (e0»  e<f>‘  Ke *  K<j)  (4.2) 

must  be  directed  along  the  outward  normal  to  the  yield  surface  at 
the  stress  point.  Thus 

q  =  X  Vf  (4.  3) 

AS** 

where  X  is  a  non-negative  scalar.  Since  X  is  unknown  Eq.  (4.  3) 
is  equivelent  to  three  scalar  equations  in  terms  of  the  unknowns 
nfl,  n^,  m0,  ny  v^,  vn,  as  required. 

At  some  points  of  the  yield  surface*  the  normal  may  not  be 
uniquely  defined.  If  the  stress  point  is  in  a  "crease"  formed  by  the 
intersection  of  two  smooth  surfaces 

f  1  *  1  f2  =  1  (4. 4) 

then  the  strain-rate  vector  can  be  any  combination  with  non -negative 
coefficients  of  the  normals  to  the  two  surfaces  forming  the  crease: 

q  =  Xj  Vfj  +  X2  Vf2  (4.  5) 

Observe  that  in  this  case  Eqs.  (4.  4)  for  the  yield  condition  provide  two 
constitutive  equations;  since  X j  and  X2  are  now  independent  unknowns* 
(4.  5)  provides  only  the  necessary  two  remaining  equations. 
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In  degenerate  cases  the  above  reasoning  may  be  generalized 
to  higher  order  intersections  which  result  in  more  yield  conditions 
of  the  form  (4. 4)  and  more  unknown  X's  in  the  flow  law.  However, 
we  observe  that  Eqs.  (3.  22)  and  (4.  4)  provide  five  equations  con¬ 
taining  only  rig  ,  n^f  and  s  as  unknowns,  and  that  the  velocity 

unknowns  v^  and  v^  occur  only  in  the  two  independent  equations  of 
(4.  5).  Therefore,  any  hypothesis  which  added  to  (4.  4)  at  the  expense 
of  (4.  5)  would  generally  lead  to  over -determined  stresses  and  under¬ 
determined  velocities  and  hence  be  unsolvable. 

The  above  feasdftlhg  has  an  Interesting  result  in  the  case  of 
a  piecewise-linear  yield  surface.  If  (4.  1)  represents  a  linear  surface, 
then  Vf  will  be  a  constant.  It  follows  that  (3.  22)  and  (4.  1)  will  provide 
four  equations  for  five  stress  variables  whereas  (4.  3)  will  provide  three 
equations  containing  only  v^  and  vn-  Therefore,  except  possibly  in 
certain  degenerate  cases,  if  the  yield  condition  is  piecewise  linear,  the 
constitutive  equations  must  be  of  the  form  (4. 4)  and  (4.  5). 


5.  PIECEWISE  LINEAR  YIELD  CONDITIONS 


Based  upon  previous  work  by  Hill  [  20]  and  Ivlev  [21], 
Haythornthwaite  [  22]  has  shown  that  any  material  yield  condition 
may  be  conveniently  bounded  by  two  piecewise  linear  material  yield 
conditions  when  the  only  physical  information  is  a  single  test  measure¬ 
ment.  We  consider  a  representation  in  principle  stress  space  as 
shown  in  Fig.  5  and  suppose  the  point  A  to  be  a  measured  tensile 
yield  stress  o^.  Following  Haythornthwaite,  we  assume  only  that  the 
material  is  isotropic  with  equal  tensile  and  compressive  yield  stresses, 
that  it  satisfies  Drucker's  stability  postulate,  [8,  9]  ,  and  that  the 
sharply  defined  rigid-perfectly  plastic  yield-point  load  represents 
useful  information.  It  follows  that  the  yield  surface  must  be  convex 
with  symmetry  every  30*.  Therefore,  if  the  tensile  yield  stress 

is  known,  the  yield  surface  must  pass  through  point  A  in  Fig.  5  and 
must  lie  between  the  solid  and  dashed  curves.  The  inner  solid  curve 
is  Tresca's  [23]  condition  of  maximum  shearing  stress 

max  [  |<Tj  -<72l  ,  |cr2-a3l,  | o\j  -o^  I  ]  =  (5.1) 

whereas  the  outer  dashed  curve  is  Hill's  [20]  condition  of  maximum 
reduced  stress 

max  lor.  -  (o^  +  cr2  +  <t3)/3|  =  <r0  1  =  1,2, 3  (5.2) 

In  the  analysis  of  thin  plates  or  shells  where  the  state  of  stress 
is  approximately  plane,  the  information  in  Fig.  5  may  be  represented 
in  a  ffj,  or2  plane.  For  flat  plates  under  bending,  the  yield  relation  be- 
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tween  principal  bending  moments  is  the  same  as  that  between  principal 
stresses(  so  that  Haythornthwaite  [22]  was  able  to  use  the  bounding 
surface  lemma  to  bound  the  yield-point  load  of  a  circular  plate  under 
uniform  pressure* 


If  both  bending  moments  and  direct  stresses  are  present,  as  is 
generally  the  case  for  curved  shellsi  the  yield  surface  will  no  longer 
have  the  simple  character  it  does  in  a  a j,  e...  space.  Rather,  four  stress 
dimensions  will  be  required,  and,  due  to  the  differences  in  integration 


generally  be  non-linear  even  if  the  material  yield  condition  is  piecewise 
linear. 


One  method  of  regaining  a  piecewise-linear  problem  is  to 
approximate  the  uniform  shell  by  an  idealized  sandwich  one.  An  idealized 
sandwich  shell  consists  of  two  sheets  each  of  thickness  t',  separated 
by  a  core  of  thickness  h*.  The  sheets  each  have  a  tensile  yield  stress 
Oq1  and  carry  no  shear;  the  core  has  no  in-plane  strength  but  is  sufficient 
ly  strong  in  transverse  shear  to  maintain  separation  of  the  sheets.  The 
sheet  thickness  t'  is  assumed  sufficiently  small  so  that  the  stresses  do 
not  vary  across  each  sheet. 

If  cr ^  (0=  1,  2)  denote  the  principal  stresses  in  the  top  and  bottom 
sheet  respectively,  the  resultant  forces  and  moments  are 

N/J  =  +  ff/3+)  *'  M/T  *(V  '  *0  >  h,t'  (5, 3) 

Evidently  Eqs.  (5.  3)  can  be  solved  for  the  four  stresses  as  linear 

functions  of  the  forces  and  moments  .  It  follows  that  any 
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material  yield  condition  which  is  piecewise  linear  in  the  stresses 

will  be  piecewise  linear  in  N.  and  M.  .  Therefore,  it  can  be  re- 

P  P 

presented  as  a  convex  polyhedron  in  a  four  dimensional  generalized 

stress  space  with  axes  N.  and  M.  . 

P  P 

On  the  other  hand,  for  a  uniform  shell  of  thickness  h  and  yield 
stress  <Tq,  the  stress  resultants  are 

h/2  h/2 

N.  =  \  <Tn  dz  M  =  \  or  zdz  (5.4) 

P  J  P  P  J  P 

-h/2  -h/2 

At  least  in  theory,  Eqs.  (5.  4)  can  be  solved  for  the  stresses  and  sub¬ 
stituted  into  the  material  yield  condition  to  obtain  the  yield  surface  in 
generalized  stress  space.  Whether  the  material  yield  condition  is 
piecewise  linear  or  not,  the  yield  surface  will  be  convex  but  will 
not  generally  be  a  polyhedron. 

We  wish  to  investigate  the  relation  between  the  yield  surface  of  the 
uniform  shell  and  the  yield  polyhedron  of  the  sandwich  shell.  To  this  end 
we  first  note  that  if  a  single  force  or  single  moment  is  acting  and 
produces  yield,  then  its  magnitude  is 

t  i  i 

Nq  =  2oq  t 

for  the  sandwich  shell  and 

>V°bh 

for  the  uniform  shell. 

The  usual  method  of  choosing  the  sandwich  shell  parameters 


i  i  i 

M0  =  Oq  k  t 


1^  =  (1/4)  <rQh^ 


No'  aI*> 


is  so  that  [  7]  , 


Mq  *^0 


(5.5) 


Now,  it  is  evident  that  the  yield  polyhedron  is  fully  specified  by  its 


vertices  and  that,  regardless  of  the  yield  condition,  each  vertex  is 

i  i 

specified  by  giving  values  of  the  four  pure  numbers  . 

Therefore,  if  (5.  5)  holds  and  if  the  same  material  yield  condition  is 
used  for  the  sandwich  and  uniform  shells,  all  vertices  of  the  yield 
polyhedron  will  lie  on  the  yield  surface.  Since  the  polyhedron  and 
yield  surface  are  both  convex,  it  follows  that  the  polyhedron  must 
lie  inside  the  yield  surface*  Therefore,  the  yield  polyhedron  for  the 
sandwich  shell  defined  by  (5.  6)  will  provide  a  lower  bound  on  the 
yield-point  load  for  a  uniform  shell  with  the  same  material  yield 
condition. 


A  different  result  is  obtained  if  we  visualize  taking  the  uniform 
shell,  leaving  its  thickness  unchanged,  but  compressing  all  of  its 
material  into  its  top  and  bottom  surfaces.  Thus 

h'  =h  o-0't'a|a0h  (4.7) 

whence 

>V  V"21^  (4-8> 


For  any  state  of  stress  this  sandwich  shell  will  be  as  strong  as  the 
uniform  shell  in  resisting  moments.  It  follows  that  the  resulting  yield 
polyhedron  must  lie  outside  of  the  yield  surface.  Therefore,  the  sand¬ 
wich  shell  defined  by  (5.  7)  will  provide  an  upper  bound  for  a  uniform 
shell  with  the  same  material  yield  condition. 
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It  is  convenient  to  define  dimensionless  generalized  stresses 
in  terms  of  the  parameters  of  the  uniform  shell  by 
nf)  =  N/S  /N0  =  +  up  )/(2<Tq  ) 

mp*  lAp  / Mq  =  (cr^~  -  ^+)/(2ko0')  (5.  9) 

where 


k  =  1  for  the  shell  of  ( 5.  6) 
k  =  j  for  the  shell  of  (5.7) 


(5.  10) 


The  stresses  can  then  be  written  in  the  form 

.  .....  ....  I?#.  . .  -  l-li  . 

Substitution  of  (5.  11)  into  the  yield  condition  will  give  a  circumscrib¬ 
ing  polyhedron  for  k  a  i  and  an  inscribing  one  for  k  =  1. 

An  alternative  circumscribing  polyhedron  can  be  obtained 
without  reference  to  a  sandwich  shell.  We  observe  that  if  the  stresses 
in  a  uniform  shell  must  satisfy  a  set  of  linear  inequalities  of  the 
form  JLi£  0 ,  then  the  generalized  stresses  must  satisfy 

Lj(n^0  L^m^JsO  (5.12) 

together  with  further  inequalities  which  represent  the  interaction  be¬ 
tween  force  and  moment.  Since  the  yield  surface  is  constrained  by 
(5.  12)  as  well  as  further  inequalities*  it  follows  that  the  convex  poly¬ 
hedron  (5.  12)  must  circumscribe  the  yield  surface.  We  shall  refer 
to  it  as  the  limited  interaction  polyhedron. 

For  plane  stress  Oj  =  0*  Tresca's  yield  condition  (5.  1)  reduces  to 
max  [  lOjI,  I  cr2l,  kj-CTj  |  ]  /<rQ'  ■'  1 


(5.  13) 
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This  must  be  satisfied  by  both  the  top  and  bottom  sheets  of  the  sandwich 
shell.  Therefore*  it  follows  from  (5.  11)  with  k  «*  1,  that 

max  [  I  n^l  +  I  m^l  ,  I  nj-n2  I  +  I  rrij  -  m2l]  =  1  (5.  14) 

is  a  lower  bound  for  any  material  yield  condition  with  tensile  yield  stress 

°o  =  V 

For  the  maximum  reduced  stress  condition  with  =  0  we  obtain 
max  [  Iff, -i^l  .  =  lj*  »/3=  1,  2;«^/3  (5.15) 

whence,  in  view  of  (5.  11)  with  k  =  | 

max  [  I  n^  -  ±n^  I  +  j  I  m^  -  im^  I  ,  \  I  nj+n2  I  +  ^  I  mj  +m£  I  ]  =  1 

(5.16) 

The  limited  interaction  polyhedron  in  this  case  is 

max  [  I  n^-jn^  |  ,  £|nj+n2l,  jlmj+m^J]  =1 

(5.  17) 

Either  of  (5.  16)  or  (5.  17)  provide  upper  bounds  for  any  material  yield 
condition. 

As  discussed  in  [  22]  and  [  24]  ,  a  similar  analysis  in  which  the 
roles  of  Tresca's  and  Hill's  material  conditions  are  essentially  inter¬ 
changed  may  be  carried  out  for  the  case  in  which  the  single  experimental 
number  is  the  shearing  yield  stress. 


6.  EXAMPLE:  CYLINDRICAL  SHELL 


For  an  axially  symmetrically  loaded  circular  cylindrical  shell 
we  choose  the  axial  direction  x  and  circumferential  direction  0 
as  principal  directions  1  and  2  respectively.  According  to  the  usual 
assumptions  of  thin  cylindrical  shell  theory*  m ^  =  m^  is  a  reaction 
and  can  be  eliminated  from  the  problem.  Further*  if  axial  load  is 
applied  only  at  the  end  of  the  shell*  the  axial  stress  n^=  t  is  constant 
along  the  shell  and  can  be  regarded  as  a  parameter.  Thus*  we  can 
reduce  the  various  polyhedra  to  polygons  in  an  n^  =  n*  m  _  =  m  space. 

We  consider  first  the  form  of  the  Tresca  polygon.  The  inequali¬ 
ties  implied  by  (5.  14)  can  be  written 


(6.  1) 


(6.  2) 


In  (6.  2),  the  value  of  m^  is  of  no  concern  and  it  is  necessary  only  to  be 
sure  that  each  left-hand  side  is  less  than  each  right-hand  side.  Eliminating 
the  redundant  and  tautological  members  of  the  resulting  set  of  16  inequal¬ 
ities  plus  the  4  inequalities  of  (4«  1)  we  obtain  the  12  inequalities 


"  }■■•(' 

-1+t  J  J  1+t 

-1+1  tl  s  m  s  1-  I  tl 

"2+t  s  2n±mS±2+t 


(6.3) 
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Inequalities  (6.  3)  define  a  polyhedron  in  an  n»  m,  t  space  which  is 
symmetric  with  respect  to  the  m  axis  and  the  plane  m  ■  0,  and 
which  inscribes  the  actual  yield  surface  when  the  tensile  yield  stress 
is  known.  The  specific  polygons  for  positive  constant  t  and  positive 
m  are 


n  =  -  1  +  t 


[2-t+2n 
1-t 

Z+t-Zn 


n  a  1 


i  ^  t  ^  1 

n  =  -  1  +  t 


m  =  1  - 1  ( 6. 4) 


n  =  1 


The  heavy  curves  in  Fig.  6show  typical  polygons  (6.4). 

A  similar  analysis  may  be  made  for  the  maximum  reduced 
stress  polyhedron.  Elimination  of  m^  from  (5.  16)  leads  to 
-8+6 1 1 1  ^  3m  <  8-61 1 1 


12  -  2t 

/ 

-8  +  2t  ^  4n±3m£  8  +  2t 


(6.  5) 
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For  an  upper  bound  on  the  tension-test  measured  shell*  we 
supplement  (6. 5)  with  the  inequalities  obtained  from  (5.  17) 
to  obtain  the  polygons 

0  st  =£2/3  2/3  St  a  1 


2n  *  -  4  +  4t 

r 

12  -  10t  +  8n 
3m  =|s  -  6t 

12  -  2t -  8n 

2n  =  4  -  2t 


(6.6) 


I  £  t  ^  4/3 
2n  =  -  4  +  4t 
3  m  =  8  -  6t 
2n  =  4-2t 

shown  by  the  light  curves  in  Fig.  6. 

The  bounding  yield  curves  for  a  cylindrical  shell  have  been  applied 
to  the  problem  of  a  cantilever  shell  under  internal  pressure  [24]  .  Figure 
7  shows  the  resulting  bounds  on  the  yield-point  pressure  as  a  function  of 
the  dimensionless  parameter 

U  =  lH  ah 

where  A  is  the  radius  and  2H  the  shell  thickness.  Similar  results  when 
the  shearing  yield  stress  is  known  are  also  given  in  [24]  . 


7.  LIMITATIONS  AND  EXTENSIONS 


The  aim  of  the  preceding  sections  has  been  to  give  the  theo¬ 
retical  background  of  a  practical  theory  of  plastic  shells.  Thus,  we 
have  attempted  neither  to  give  a  complete  catalogue  of  available 
problem  solutions,  nor  to  give  the  most  general  possible  theory. 
With  regard  to  problem  solutions,  a  representative  selection  and 
extensive  bibliography  may  be  found  in  [  10]  . 


Among  the  physically  present  physical  concepts  which  have 


elastic  and  plastic  strains,  the  effect  of  strain -hardening,  changes  in 
stress  distributions  due  to  small  geometry  changes  induced  by  the 
loads,  and  nonlinearities  due  to  finite  strains.  Hodge  and  his  associates 
[  25,  26,  27,  28]  have  considered  various  simple  problems  in  which 
elastic  strains  and  strain-hardening  were  included.  In  every  case 
investigated,  a  representative  load -deformation  curve  had  the  qualita¬ 
tive  form  shown  in  Fig.  8  [25]  .  Based  upon  these  examples,  it  appears 
reasonable  to  assume  that  structures  made  of  reed  material  exhibit 


qualitatively  different  behavior  depending  upon  whether  the  load  is  above 
or  below  the  yield-point  load  of  the  same  structure  made  of  an  idealized 
rigid/ perfectly-plastic  material.  Thus,  if  the  purpose  of  the  investigation 
is  to  determine  only  the  load  value  at  which  this  qualitative  difference 
occurs,  the  rigid/perfectly-plastic  model  considered  herein  will  provide 
a  reasonable  estimate  for  the  desired  information.  Further,  Figure  8 
indicates  that  if  more  detailed  information  is  desired  for  loads  less  than 


the  yield-point  load,  it  is  reasonable  to  neglect  strain-hardening  and 
consider  an  elastic/perfectly-plastic  material.  Finally,  for  information 
at  loads  above  the  yield-point  load,  elastic  strains  are  relatively  un¬ 
important  and  one  may  use  a  rigid/ strain -hardening  material  as  a  model. 

The  situation  with  regard  to  small  and  large  geometry  changes 
is  less  clear.  It  has  been  shown  by  Haythornthwaite  [  29]  that  the  be¬ 
havior  of  a  beam  whose  ends  are  fully  fixed  as  to  both  slope  and 
separation  is  quite  different  from  one  whose  ends  are  clamped  but 
free  to  move  towards  each  other  under  load.  For  the  former, 
even  a  small  deformation  of  the  order  of  half  the  beam  height  intro¬ 
duces  axial  forces  which  substantially  raise  the  load-carrying  capacity 
of  the  beam.  Similar  results  for  circular  plates  were  found  by  Onat 
and  Haythornthwaite  [  30]  . 

The  beam  and  circular  plate  problem  have  proved  solvable 
because  of  the  fact  that  they  deformed  into  easily  characterized  simple 
shapes.  For  other  shell  problems,  the  initial  velocity  field  at  the 
yield-point  load  predicts  that  elementary  shell  shapes  such  as  spheres, 
cones,  or  cylinders,  deform  into  complex  shapes  which  are  not  easily 
characterized.  Therefore,  a  general  theory  of  the  post-yield  behavior 
of  shells  must  probably  await  a  more  general  technique  for  determin¬ 
ing  the  yield -point  load.  In  view  of  the  complexity  of  solutions  for  such 
simple  shapes  as  spheres  or  cones,  it  appears  almost  certain  that  any 
general  approach  must  be  primarily  numerical.  A  first  step  in  this 
direction  has  been  taken  by  Onat  and  Lance  [  16]  ,  for  a  shallow  conical 
shell,  but  it  is  not  yet  clear  if  their  methods  can  be  generalized. 
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With  regard  to  truly  large  deformations*  the  small-strain 
theory  presented  herein  is  wholly  inadequate*  However*  for  very 
thin  shells*  it  does  appear  reasonable  in  this  case  to  neglect  bending 
stresses  entirely  and  construct  a  membrane  theory  of  finite  shell 
deformation.  This  has  been  done  by  Salmon  [31]  for  an  initially 
cylindrical  shell.  Obviously  more  work  remains  to  be  done  in  this 
area  also. 
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Figure  2.  Yield  curve  for  rectangular  beam 


Figure  7.  Bounds  on  yield-point  load  for  cylindrical  shell 
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